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Abstract 

The Cervera et al. formula, the best known approximate formula of neutrino oscillation proba- 
bility for long-baseline experiments, can be regarded as a second-order perturbative formula with 
small expansion parameter e = Amli/Amf^ — 0.03 under the assumption S13 ~ e. If #13 is large, 
as suggested by a candidate v e event at T2K as well as the recent global analyses, higher order 
corrections of S13 to the formula would be needed for better accuracy. We compute the corrections 
systematically by formulating a perturbative framework by taking #13 as S13 ~ yfi ~ 0.18, which 
guarantees its validity in a wide range of #13 below the Chooz limit. We show on general ground 
that the correction terms must be of order e 2 . Yet, they nicely fill the mismatch between the ap- 
proximate and the exact formulas at low energies and relatively long baselines. General theorems 
are derived which serve for better understanding of 5-dependence of the oscillation probability. 
Some interesting implications of the large #13 hypothesis are discussed. 
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I. INTRODUCTION 



One of the most important progresses in particle physics in the last decades is the discov- 
ery of neutrino masses [l| and the lepton flavor mixing 0]. It was done through observing 
neutrino oscillation phenomena and it constitutes, up until this moment, only available ex- 
perimental method for measuring lepton mixing parameters. Am^ and #23 are determined 
by atmospheric neutrino observation by Super-Kamiokande and then by accelerator 

neutrino experiments 0, 0]. Am^ and 6*12 are measured independently by two types of 
experiments, the KamLAND reactor experiment [8-10] and the solar neutrino observation 
using various experimental techniques. For the latest results and for a review of the solar 
neutrino experiments see e.g., 0,0 and 0, respectively. The remaining mixing angle 9 



13 



is being explored by the ongoing and the upcoming accelerator [14], [15] and reactor neutrino 



experiments [3-18]. If it turned out that 13 is not too small, we may proceed to measure 
CP violation by the lepton Kobayashi-Maskawa (KM) [l9| phase 84, to which we refer just 
8 in this paper. 

It is expected that precision measurement is required to determine 5 because CP violation 
effect is tiny due to suppression by the two small factors, AmS^/Am^ and the Jarlskog 



coefficient J = C12S12C23S23C13S13 2(f . Therefore, understanding of full complexity of neutrino 
oscillation phenomena would be of some help e.g., to design future experiments. An example 
of such is the parameter degeneracy 2lM23l| , the problem of multiple copy of the solutions 
of mixing parameters allowed by given sufficient but limited numbers of experimental data. 
See [24] for a comprehensive overview of this phenomenon. To facilitate understanding of 
qualitative features of the neutrino oscillation, it is crucially important to have analytic 
formula, albeit approximate, for the oscillation probability. For relatively short baseline 
experiments, such as low-energy superbeam (25-13], the matter perturbation theory works 



28l . [29j. So far, most of the analyses for long baseline of L > 1000 km were done by using 
the well known Cervera et al. formula [30]. 1 

A simple way of deriving the Cervera et al. formula is to expand the exact oscillation 
probability by small expansion parameters, e = Amli/AmJi and S13 = sin #13, both to 
second order. While the former is known to be small e ~ 0.03, the latter can be larger 
by almost an order of magnitude; Currently, it is only bounded from above by the Chooz 



limit, Si 3 < 0.18 [34M371] . We note that there is an indication for nonzero #i 3 comparable 
to the Chooz limit from global fit of the solar, reactor, atmospheric, and the accelerator 
experiments Though the statistical significance of the indication is not high 

enough to be compelling, it certainly gives a good motivation for examining effects of such 
large values of #13 . Recently, a candidate event for v e appearance has been seen in the T2K 
experiment [43[ which further strengthens the motivation for taking the large-#i 3 hypothesis 
seriously. 

If #13 is large, higher order terms of Si 3 would be needed to achieve better agreement with 
the exact oscillation probability In this paper, we compute such higher order corrections 
of S13. To facilitate systematic computation we formulate a perturbative framework by 
assuming s 13 as large as ~ y/e, which is comparable to the Chooz limit, and by taking the 
case of long enough baseline where the matter effect is comparable to the vacuum one. We 



However, it was shown that analysis of the parameter degeneracy with the matter-perturbative formula 
has been proved to give a transparent view of the phenomenon [3ll433j , such as the decoupling between 
the degeneracies. 



2 



call this framework as the u y/e perturbation theory" as opposed to the e perturbation theory 
for the Cervera et al. formula as named in j33|. We derive the second order formulas for 
the oscillation probabilities in all channels. By taking the ansatz, si 3 ~ y/e, our formulas 
will be valid for a wide range of 6*13, e < S13 < y/e. The y/e perturbation theory has been 
formulated earlier for relatively short baseline setting [44] • 

Some characteristic features of the computed formulas prompted us to think about general 
property of the oscillation probability, which resulted in the two general theorems described 
in Sec. [TTJ For example, we show that the 5 dependent terms in the ^-related oscillation 
probabilities exist only in odd terms of S13, and conversely, all the terms odd in Si 3 have 
S dependence. Thanks to the theorems, quite (un)interestingly, we can prove on general 
ground that the correction terms to the Cervera et al. formula have to be second order in 
e, eliminating the possibility of having large corrections. Yet, we will observe that at super 
long baseline, L = 4000 km, our correction terms nicely fill a sizable gap between the exact 
oscillation probability and the Cervera et al. formula, which exists in a limited range of 
energy for # 13 comparable to the Chooz limit. 

The possible large value of #13 could generate some interesting effects. For example, 
there arise terms of order e 5//2 consisting solely of 5-dependent terms in the v e appearance 
oscillation probability, which is smaller only by a factor of ~ 5 compared to the existing 
terms in the Cervera et al. formula. Another intriguing feature arises when the non-standard 
interactions (NSI) of neutrinos are included into the system. For an extensive list of the 
references for NSI including the original ones see, e.g., [45]. Some of the NSI dependent 
terms get enhanced by large #13, and decoupling of some NSI elements from the z/ e -related 
appearance probabilities no more hold. On the other hand, the smallness of the correction 
terms to the Cervera et al. formula has a consequence that, quite naturally, their influence 
to parameter determination including the issue of parameter degeneracy is quite limited. 

Following the discussion of general property of oscillation probability in Sec. [Til we for- 
mulate our y/e perturbation theory (Sec. IIIII) . and derive the expressions of the oscillation 
probabilities valid to order e 2 in v e — > and — > channels (Sec. HVj) . The charac- 
teristic features of the formulas are shown to be understood thanks to the theorems given 
in Sec. HU The accuracy of our formula is checked against the exact results for a large 
value of 6 I 13 (Sec. |V]). Then, we calculate order e 5//2 terms of the oscillation probability in 
the v e — > Up channel which could be relevant for measurement of 5 in its extreme precision 
(Sec. IVip . Implications of smallness of the large #13 corrections are explored by treating 
the parameter degeneracy (Sec. IVIIj) . The second-order formula of the oscillation proba- 
bility in the ^-related appearance channels is derived for systems with NSI (Sec. IVIIIj) . 
Concluding remarks are given in Sec. IIXI Appendices are devoted for a proof of a general 
theorem for suppression of CP phase effect (Appendix [A]), the explicit expressions of the S 
matrix elements in the #23- r otated basis (Appendix [B]) , and explicit forms of the oscillation 
probabilities in the remaining channels (Appendix IC]) . 



II. GENERAL FEATURES OF NEUTRINO OSCILLATIONS IN MATTER 

Before constructing perturbative framework of neutrino oscillation, let us discuss some 
generic features of the oscillation probability. They are interesting by themselves, and state- 
ments about cos 5 dependence do not appear to be explicitly spelled out in the literature to 
our knowledge. It will also help us to understand characteristic features of the perturbation 
theory of neutrino oscillation to be discussed in the rest of this paper. An example of this 
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is that CP phase effect appear only in terms of half-integral order of the small expansion 
parameter e in our yfk perturbation theory. 

We note that the general features of neutrino oscillation probability which exactly hold 



with |46|, |42| or without |48| constant matter density approximation can be used to prove 
some useful properties of the oscillation probabilities. That is, we call the readers' attention 
to the following two theorems: 

• Theorem A: In z/ e -related oscillation probabilities the 5- dependence exists only in odd 
terms in s 13 . Conversely, all the odd terms in s 13 are accompanied with either cos 5 or 
sin 5. 

• Theorem B: ^-dependent terms in the oscillation probability in matter, not only sin 5 
but also cos 5 terms, must come with the two suppression factors, . and the reduced 

Am 31 

(or "not quite") Jarlskog coefficient J r = C12S12C23S23S13. 

In talking about S13- and (^-dependences we assume the standard form of the lepton flavor 
mixing matrix, the MNS matrix [2| given in ([7]), the convention which will be used throughout 
this paper. 

To prove the theorem A we note that Si 3 and 5 enter into the Hamiltonian through the 
single variable z = Si 3 e l<5 . Therefore, the oscillation probability P can be written as a power 
series expansion as P = Y^ 3 f n m z n (z*) m , where f nm = f^ n for reality of P. On the other 
hand, the result obtained in 1481 by extending discussions in (49[ says that there are only 



cos S and sin 5 terms in u e related oscillation probabilities, and no higher harmonics of S. 
It means that only the terms that satisfy m = n ± 1 survive. It leaves the unique form of 
the oscillation probability P = K(s1 3 )si3 cos 5 + M(s^ 3 )si3 sin 5, where K and M are some 
functions. At the same time, this construction guarantees that all the odd terms in s 13 are 
accompanied by 5. This is nothing but the theorem A. 

Now, let us discuss the theorem B. If the adiabaticity holds existence of the suppression 
factors by the Am 2 ratio and the angle factors for sin S terms is well known. In fact, one 
can show that the latter is precisely the Jarlskog coefficient J = C12S12C23S23C13S13 thanks to 



the Naumov identity [5Cj, 151 



Am^AmlgAm^ J sin 8 = Am^Am^Am^ J sin 5, (1) 

where the quantities with over-bar in the left-hand side denote the ones in matter which 
correspond to the one in vacuum in the right-hand side of ([I]). Notice that cos 5 does not 
contain sin<5 j47[, and the proportionality between sin 5 and sin 8 is also guaranteed by the 
Toshev identity c 2 3S 2 3sin5 = c 2 3S 2 3sin(5 [52|. The equation ([TJ indicates that the sin 5 term 
vanishes in the limit of one of Am| — > 0, or vanishing limit of one of the mixing angles. 

Whether the same statement apply to cos 5 terms or not, or if the theorem B is valid 
at all for cases with generic matter density profiles, is not obvious. Apparently no general 
statement has been made in the literature. Vanishing of any 5 dependence (including cosine) 
in the absence of Am 21 can be proved similarly as the phase reduction theorem given in [33| 
as a special case of turning off the non-standard interactions. For constant matter density 
proportionality of cos S terms to the Jarlskog coefficient J (reduced coefficient J r ) in the 
z/ e -related channels (oscillation channels in the — v T sector) is explicitly proved in [47] by 
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deriving the exact forms of the oscillation probabilities. 2 Therefore, what is left is to show 
that the last statement holds under arbitrary matter density profile without recourse to the 
assumption of adiabaticity. The proof for this general case is described in Appendix [A] 



III. FORMULATING LARGE-6»i 3 PERTURBATION THEORY 

In this section, we formulate a large-(9 13 perturbation theory of neutrino oscillation. We 
use an ansatz 

/\td.^ 

s 13 ~ y/i, e= ~ 0.03 (2) 

to formulate our perturbative framework, hence the name u \/e perturbation theory". It 
implies sin 2 2# 13 ~ 0.12, the value comparable with the Chooz limit. We work in antic- 
ipation of long baselines of several thousand kilometers, so that = a/Am^ ~ 
and Am\ x L/E ~ 0(1) to be not far from the first oscillation maximum. Hereafter, 
a = 2\f2GpN e {x)E is a coefficient for measuring the matter effect on neutrinos propa- 
gating in medium of electron number density N e (x) [53], where Gp is the Fermi constant 
and E is the neutrino energy. In perturbative calculation of the oscillation probability we 
take constant electron number density approximation. 

If we take the ansatz s 13 ~ e (which corresponds to sin 2 2#i 3 ~ 4e 2 ~4x 10~ 3 ) instead 



of (|2J), we obtain the widely used Cervera et al. formula [30j by keeping terms to order 
e 2 . Of course, there may exist the other small (or large) parameters, such as Am 3l L/E 
and a/Amj! at short baselines, depending upon the experimental settings. The similar y/e 
perturbation theory for shorter baselines is discussed by assuming ta ~ \fe [44j |. 

The remaining potentially small parameters would be 7r/4 — ^ 23 , but we do not take it 
as an expansion parameter for two reasons: 3 (1) A rather large range is currently allowed 
for 623, and moreover the situation will not be changed even with the next generation 
experiments j57[. (2) As will become evident in formulating our perturbative framework # 23 
is an "external parameter" which is irrelevant in doing perturbative computation. 

We follow 33|| to formulate the perturbative treatment of neutrino oscillation. The S 



matrix describes possible flavor changes after traversing a distance L, 

y a (L) = Sapvpify, (3) 

and the oscillation probability is given by 

P{up^u a -L) = \S aP \ 2 . (4) 

When the neutrino evolution is governed by the Schrodinger equation, i4-v = Hi/, S matrix 
is given as 

■L 



S = Texp 



i I dxH(x) 



o 



(5) 



2 The difference of whether the suppression factor is provided by J or J r does not make difference in our 
following discussions. Hence, we do not discuss this point further apart from making a brief comment in 
Appendix |XJ 

3 Nonetheless, one can of course further expand the oscillation probability in terms of 7r/4 — #23 assuming 



it small, as is done e.g., in [541456 1 
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where T symbol indicates the "time ordering" (in fact "space ordering" here). The right- 
hand side of f$5§ may be written as e~ lHL for the case of constant matter density. In the 
standard three-flavor neutrinos, Hamiltonian is given (with a = 2\/2G F N e E) by 



H 















a(x) 





o" 







Am 2 21 
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where Am 2 - 
notation as 



m 2 — m 2 . In 

j 1 



(6) 



U is the MNS matrix and can be written in the standard 
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(7) 



with the notation = sin^- etc. and 5 being the lepton KM phase. 

To formulate perturbative treatment it is convenient to work with the tilde basis defined 
as v a = (Ul 3 ) a pUp, in which the Hamiltonian is related to the flavor basis one as 58 



H 



U\ 3 HU 23 



(8) 



where U 23 is defined in (J7J). The S matrix in the flavor basis is related to the S matrix in 
the tilde basis as 



where S(L) = Texp —i dxH(x) 
is given by 



S(L) = U 23 S{L)U\ 3 (9) 
The unperturbed part of the tilde-basis Hamiltonian 



H — A 



r A 

1 



(10) 



where A = A ^ 31 and r A 



2E 



Am: 



9l -t. While the perturbed part is written with another simplified 



notation r A 



"a — t~ as 

Am;, 



#1 



A 



- A 



s 13 e- iS 

s l3 e iS 






+ A 



r A S? 2 + s 13 r ACl 2 Si2 



rACl 2 Si2 








rAC 2 12 
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— s 



13 



Ar A 



{rAS 2 12 + \s\ 3 ) s 13 e iS r A ci 2 si 2 si 3 e* 

S 12 S 13 | c 12Sl2Si 3 

\c 12 s 12 s 2 n 
-s 2 2 s 2 3 



(rAS 2 2 + ^sj 3 ) s 13 e 
r A c 12 s 12 s 13 e~ lS 




i5 



(11) 



The first, second, third, and the fourth terms in ffTTj) are of order ea, e 1 , ea, and e 2 , respec- 
tively. 
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To calculate S(L) perturbatively we define Q(x) as fl(x) = e lH ° x S(x), which obeys the 
evolution equation 



i—n(x) = HMx) 
ax 

where 

Hi = e iH xjj ie -iH x 

Then, Q(x) can be computed perturbatively as 



(12) 



(13) 



tt(x) 



1 + 



" u i(x") 



+ 



rx rx px 

-i) I dx' 'H x (x') + {-if j dx'H^x') I dx"H 
Jo Jo Jo 

[ dx'H^x') [ X dx"Hi{x") [ X dx'" H^x 1 ") + 0{ 
Jo Jo Jo 



(14) 



where the "space-ordered" form in (iMj) is essential because of the non-commutativity be- 
tween H\ of different locations. Having obtained Q(x) S matrix can be written as 



S(x) 



D -iH x 



fl(x) 



(15) 



The results of S(x) matrix elements to second-order in e are given in Appendix [B] Then, 
the S matrix can be computed by making a rotation in (2-3) space S — U23SU23 ^ ^ (191) ■ 
(See ( 1B9I) .) Finally, the oscillation probability can readily be obtained by using (@J. For 
example, the one in the v e — > channel can be given by using the S matrix elements as 



P(u e = \S, e \ 2 = sl 3 \Sil /2) (-t)\ 2 + 2c 23 s 23 Re [§§M (S) 

+ 4l^)(-5)| 2 + 2^ 3 Re Is^i-sys^i-s) 



(16) 



IV. PERTURB ATIVE EXPRESSION OF THE OSCILLATION PROBABILITY 

In this section, we present perturbative expressions of the oscillation probabilities in the 
v e related and the — v T sectors to order e 2 in our y/e perturbation theory All the formulas 
in this section are given in the form as 

in which we use L to denote the baseline distance. Some order e 5//2 terms will be discussed 
in Sec. I VII The antineutrino probability can be obtained from the neutrino probability by 
the replacement as P{p a —> V(i\ 8, a) = P{y a vp] —5, —a). Similarly, the T-conjugate one 
is given by P{y^ —> v a \ 6, a) = P{y a — > —5, a). Remember that the following abbreviated 

notations are used: A = A !^ 1 , r& = , and = . % . 
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A. Oscillation Probabilities in v e related sector 



There is no zeroth oder term in the appearance channels. With use of the reduced Jarlskog 
coefficient J r = C12S12C23.S23.S13, the order e 1 , e 3//2 , and e 2 terms in the oscillation probability 
in the v e — > channel are given by 

- r A ) AL 



p(3/2) 



P(l) 



As 2 s 2 - 



l-r A ) 



sin 



rA 



r A (l - r A ) 



cos ( 5 — 
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. r A AL . (1 — r A )AL 
sm sin . 



(18) 
(19) 



p(2) 
en 
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4C 23 C 12 S 12 



Is 



23 



l+r A ) 2 



sm 



r A AL 



'13 



1 



+ 2s 



23 



2s 



■ r A ) 
r A 



13" 



2s 2 s 2 - 



2 2 
S 12 S 13" 



7"A?" A 

l-r A ) 
'a 



sm 



;i - r A )AL 



(AL) sin(l — r A )AL. 



(20) 



l-r A ) 3 " id (l-r A ) 2 . 

As is obvious from flB9[) . the similar expressions Pjr for the z/ e — > u T channel can be 
obtained from pj$ in (TT8|) - (T20|) by the transformation C23 — > — S23 and S23 — >■ C23. The 
explicit expressions are given in Appendix Given pJJ and P®, P® can be obtained by 
using the perturbative unitarity relation 4 

PW = (5- n - P W - P W f 



0,1,3/2,2). 



Therefore, we do not present their explicit forms. Notice that Pi^ 2 "* + Pir^ 



(21) 
as it should 



because there must be no 5 dependent terms in P ee [49|, |60 



B. Oscillation Probabilities in i/„ — u T sector 



The order-by-order perturbative formulas of the oscillation probabilities in the — > 
channel to order e 2 can be computed via a similar manner. The results are given by 



p(0) 



1 - 4c 23 S 23 Sin2 



AL\ 



(22) 



PW = -4sLs 2 



, 23 c, 13 



l-r A 



sm 



1 - r A ) AL 



- 2c 2 s 2 

zc 23 6 23 



r A 



'13 



4- Ac 2 s 2 s 2 - 



1 - r A 
1 

1 - r^) 5 



( AL) sin AL 



. (1 + r A )AL . (1 — r A )AL 
sm sm . 



(23) 



4 In fact, we computed all the Pea (a — e, /i, t) by the same procedure and explicitly verified the pertur- 
bative unitarity relation. 
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P<f) = -8J r oosi(4-4,)f^-sm 2 (^) 



+ 8 J r cos 5 — - — sin 2 

r A (l - r A ) 

ir/ x 2 . AL . r A AL (1 - r A ) AL 

— lo J r cos os 23 — 7- -sin — -sin — - — cos , (24) 

r A [l — r a) 2 2 2 
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2coqS 



2 2 
23*23 
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( AL) sin Ax 
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The similar order e 1 , e 3//2 , and e 2 terms P^r in the oscillation probability in the 



(25) 



— >■ I/r 



channel are given in Appendix For ^ — >• i/ e channel, P$>(6) (i = 0,1,3/2,2) can be 

obtained from its T-conjugate as Pe l J(—5). It is then straightforward to verify that the 
similar perturbative unitarity relation, with the indices e being replaced by /1 in ( 12~TT) . holds, 
the task we have explicitly executed. 



C. Understanding some characteristic features of the perturbative formulas 

We first note that the S dependence appears only in terms with half-integral order in 
e, and all the terms of half integral order in e contains 5 dependence. It can be readily 
understood by the theorem A in Sec. |TT| because odd terms of S13 have to be half-integral 
order of e in our yfe perturbation theory. 

Despite that the Cervera et al. formula is the second order formula for small Si 3 ~ e, 
practically it is often used even for relatively large #13, for example, in the analysis of 



parameter degeneracy [24], |59(. Therefore, it is a legitimate question to ask how large the 
corrections terms to the formula can be for large #13. We therefore classify each term in 
the oscillation probabilities into the two categories, the one which exist in the Cervera et al. 
formula (denoted for brevity as the Cervera terms), and the ones which do not (non-Cervera 
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terms). The non-Cervera terms in each channel are the terms with factors of either sf 3 or 
vas± 3 . In the v e — > (u e — > v T ) channel they are the last two lines in (|20|) ( (1C3j) ). which 
come from the last term in (|T6|) . 

One notices that the non-Cervera terms arise only from the order e 2 terms. To understand 
this feature let us first note that there can be no non-Cervera terms of order e 1 in our y/e 
perturbation theory. While the e 1 suppression is provided by the factor of either r& or 
s 2 3 , they can be among the e 2 terms in the e perturbation theory, which is included in the 
Cervera terms. Therefore, the largest possible non-Cervera terms may be contained in the 
order e 3//2 terms. However, they do not exist for the following reason: The theorem A proved 
in Sec. [Til states that all the odd terms in s 13 must be accompanied either by cos 5 or sin 5. 
Then, the theorem B dictates that all the 5 dependent terms must have extra suppression 
factor ta = Am^/ Am^ ~ e. Therefore, sf 3 terms in the non-Cervera part of the oscillation 
probability are actually of order e 5//2 . As a consequence, the non-Cervera terms exist in our 
second-order formulas only with the suppression factors of s^ 3 or tas1 3 . Hence, they are of 
order e 2 , excluding the possibility of yielding large corrections of order e 3 / 2 or lower from 
the non-Cervera terms. 

V. ACCURACY OF THE APPROXIMATE FORMULA FOR LARGE 8 13 

In this section, we examine numerical accuracy of our perturbative formulas of the oscil- 
lation probabilities computed to order e 2 for #13 of the order of the Chooz limit. We focus on 
the v e — > Vfj, channel whereas we also did the similar analysis in the — > channel [61]. To 
display accuracy of the approximate formula, two typical ways of plot are available, namely, 
the absolute difference |-P e e * act — -P e 2 ^ d |, and the relative difference |i^ act — -P e 2 ^ d |/-P e e ^ act , where 
pexact an( j p2nd d erio t e; respectively, numerically evaluated exact oscillation probability and 
our approximate formula to order e 2 . 

In Figs. [TJ and [2J the absolute and the relative differences, respectively, are presented by 
using the color graduation plots in E — L space. We note that the comparison between 
the Cervera et al. formula and the exact result was performed in [58[. We use the similar 
format to make easier the comparison between our and their results. A comparison between 
the relevant panels in Figs. 3 and 4 of |58| and Figs. [1] and [2] indicates that our large-^13 
correction terms improves the accuracy of the approximate formula in a wide region in E — L 
space. It is even more so considering that our assumed value of sin # 13 is about a factor of 
two larger than their largest value. The features of the plots with other values of S are quite 
similar to those presented in Figs. [TJ and [21 and hence we do not present them. For the same 
reason, only the case of normal mass hierarchy is shown in Fig. [2j 

To see more clearly how good (or bad) are the approximations by our and the Cervera 
et al. formula, we present in Fig. [3] the oscillation probability as a function of neutrino 
energy calculated numerically (denoted as "exact", green dashed line), computed by using 
the Cervera et al. formula (blue dash-dotted line), computed with our perturbative formula 
(red solid line). The left and the right panels in Fig. [3] are for baselines L = 1000 km and 
L = 4000 km, for which the matter density is taken as 2.8 g/cm 3 and 3.6 g/cm 3 , respectively. 
The mass hierarchy is taken as the normal one, Am 2 x > 0. #13 is taken as sin# 13 = 0.18 
which is close to the Chooz limit, and S = 0. The values of the remaining mixing parameters 
used are given in the caption of Fig. [TJ 

As seen in Fig. [31 the difference between the exact and the Cervera et al. formula is 
very visible (modest) at L — 4000 km (L = 1000 km). It is notable that the higher 
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FIG. 1: The absolute difference between the numerically evaluated exact oscillation probability 
and our approximate formula to order e 2 in the v e — > channel, |-Pe^ act — P^a \i is presented 
by using the color graduation plot in E — L space. The top (bottom) two panels are for the 
normal (inverted) mass hierarchy, the left neutrino and the right anti-neutrino channels. The 
correspondence between colors and the probability difference is given at the top of the figure. The 
matter density is taken as 3.0 g/cm 3 for all baselines. #13 is taken as sin #13 = 0.18 and 5 = 0. 
The remaining mixing parameters are chosen as Am^ = 2.4 x 10~ 3 eV 2 , Am^ = 7.7 x 10~ 5 eV 2 , 
sin 2 623 = 0.5, and tan 2 #12 = 0.44. 



I I <0.001 0.001-0.01 0.01-0.05 ■ 0.05-0.1 ■ >0.1 

v e ^ V(J Ami, > 5=0 v e ^ Amf, > 5=0 




E[GeV] E[GeV] 

FIG. 2: The relative difference |-Pe^ act — P^ d I /-f^ act is presented with the same format and by 
using the same values of the parameters as in Fig. (TJ 
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FIG. 3: Comparison between the exact oscillation probability P{v e —> computed numerically 
as a function of energy (green dashed line), the one calculated by the Cervera et al. formula (blue 
dash-dotted line) , and with our formula with large #13 corrections (red solid line) . The left and the 
right panels are for baselines L = 1000 km and L = 4000 km for which the matter density is taken 
as 2.8 g/cm 3 and 3.6 g/cm 3 , respectively. #13 is taken as sin 6*13 = 0.18 and 5 = 0. The values of 
the remaining mixing parameters are the same as given in the caption of Fig. [TJ 



order corrections of S13 incorporated in our formula nicely fill the gap between them. To 
understand the nature of the gap we have examined the other cases of mass hierarchy and 
(anti-)neutrino channel. The features of the three curves presented in Fig. [3] are very similar 
to the ones of antineutrino channel with the inverted hierarchy, both = . a - 2 > 0. The 
features in the other two cases, antineutrino channel with the normal hierarchy and neutrino 
channel with the inverted hierarchy, both < 0, are similar to each other. They are not 
presented because the differences among the three curves are much smaller for both L = 1000 
km and 4000 km. We have checked that the features mentioned above are very similar for 
other values of 5. 

As was proved in Sec. IIVCI the higher order corrections (the non-Cervera terms) are of 
order e 2 ~ 10~ 3 , and hence the large gap between the exact and the Cervera et al. formula, 
in particular the one at L = 4000 km, must arise as a result of some enhancement. It is 
caused by the factor 1/(1 — r^), as seen from the relevant energies of the gap and the features 
of the other channels and hierarchies; It is nothing but the enhancement due to the MSW 
resonance 



63|. It is interesting that our formula, though perturbative, can incorporate 
the enhancement effect for baselines up to several thousand km for which AL ~ 5, not too 
far from ~ 0(1). For baseline of order L ~ 10000 km the resonance enhancement becomes 
more significant and the difference between the exact and our formula (as well as the Cervera 
et aVs) blows up. Of course, it is outside the region of validity of the perturbative treatment, 



and one need to sum up 1/(1 — ta) effect. It is beyond the scope of this paper. 



VI. P{y e -> z^) TO ORDER e 5 / 2 AND ITS 5 DEPENDENCE 

For large #13 comparable to the Chooz limit, it may be meaningful to analyze the order e 5//2 
terms, though it has an extra suppression of ^ft ~ 0.2 compared to the order O (e 2 ) terms. 
It could be particularly relevant for the u~ — > channel for which an extreme accuracy 
may be reached e.g., by neutrino factory [64]. Our interest in the e 5//2 terms is primarily 
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due to that they consist only of 5- dependent terms, the property enforced by the theorem 
A. Fortunately, it is easy to compute the order e 5//2 terms in the oscillation probability by 
using the S matrix elements to order e 2 listed in Appendix iBl 

The results of the order e 5 / 2 terms in P{y e — > v^) which is to be added to the e 1 , e 3 / 2 , 
and e 2 terms given in Sec. IIVAI reads: 



d(5/2) o t 2 r A r A , . a (1 - r A )AL 
PW ; = 8J r sf 3 — — — cos g sin 



+ 8J, 



-2s 2 rA + f c 2 - s 2 )— + s 2 rArA 

13/1 „ \2 ^ V C 12 A 12/„ ^ 6 12i 



rx(l-r A ) L 13 (1-^a) 2 

x cos I 5 — 




\ 2 J 2 

QT 2 r A ( /\ t \ fx AL \ ■ r ^ AL (1~Ta)AL 
+ 8Jrg i3 ( 1 _ rA ) 2 ( AL ) cos ( 5 - — ) Sm _ ~ C ° S 2 

rl ,. n / r A AL \ . r A AL 

■ AL cos sm 

i-r A y 1 v 2 ; 2 

r\ ( (l + rAAL\ ( 



A) \ / 

I / An / r (l+r^ALA . (l-r A )AL 
r A (l-r A ) V 2 / 

i t r A / 2 r A 2 \ / (l-r A )AL \ , (1 — r A )AL 

~ AJr r A (l-r A ) [ S ^— A ~ Sl2rA J (AL) C ° S » 1 Sm " 



(26) 



Pir^ can be obtained from PefJ 2 ^ by the transformation c 2 3 — > — s 2 3 and s 23 — > c 23 . 

The distinctive feature of order e 5//2 terms in ( 12 6 p . as compared to e 3 / 2 terms in (I19p . is 
that there exist much more profound type of 5 dependence, not only correlated (as argument 
of cosine function) with 4p but also to Va £ lL and ( 1±r ^) AL - jf measurement is sufficiently 
accurate to resolve such correlations different from the vacuum type, it must merit to achieve 
higher sensitivity to detect CP violation and accurately measure 5. However, quantitative 
analysis to examine this feature to reveal the required experimental conditions, such as 
which energy resolution and how many baselines are required etc. is beyond the scope of 
this paper. 



VII. PARAMETER DETERMINATION AND DEGENERACY 



The analysis of parameter degeneracy e.g., in |24| can, in principle, be repeated with 
our formula with higher-order corrections of S13. However, since the probability contains 
quartic terms of S13, one has to deal with eighth-order equations of S13 to obtain the full 
degeneracy solutions, a formidable task to carry out. Fortunately, we know empirically that 
apparently there is no other solution beside the known eight-fold degeneracy for # 13 below 



the Chooz limit 65] , though it may worth further examination. Therefore, in this paper we 
limit ourselves to the known degeneracy solutions and estimate corrections to them due to 
the large-#i3 correction terms. 

We first discuss relationship between the determined mixing parameters with and without 
the non-Cervera terms. Given the "observable", a pair of the oscillation probabilities P efl = 
P(u e — » Ufj) and P efl = P{p e — >• Vp) at certain neutrino energy E, we consider the problem 
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of determining the set of parameters (s = Si 3 ,<5). If we use the Cervera et al. formula, the 
observable are related to mixing parameters (so, 5c) as 

P eM - Z = X±s c + Y±s c (cos 5 C cos A 3i ± sin S c sin A 3 i) , 

P eM - Z = X T s 2 c - Y T s c (cos 5 C cos A 3 i =f sin 5 C sin A 31 ) , (27) 

where A 31 = \AL/2\, and X, Y, and Z are the coefficients whose explicit forms may be 
constructed from the formulas given in Sec. IIVI (or, see equation (2.7) in j24[). The ± signs 
in this section imply the normal and the inverted mass hierarchies. 

When the non-Cervera corrections are taken into account the same observable yield 
slightly different set (st, 5t) of mixing parameters as 

P efl — Z = X±s T + Y±st (cos 5t cos A 3i ± sin 5t sin A 3 i) + Pnc, 

P efX — Z = X T s T — Y t st (cos 5t cos A 31 =f sin 5t sin A 31 ) + Pnc, (28) 

where Pnc an d Pnc denote the non-Cervera corrections in neutrino and anti-neutrino chan- 
nels, respectively. The point is that the arguments in Pnc and Pnc is in principle (st, St) 
but it can be replaced by (so, 5c), the known quantities, because all the non-Cervera cor- 
rections are second order in e. It is then easy to compute the difference between st and sc 
and 5's to leading orders in e as 

t F T sin (5 C ± Agi) P NC + Y± sin (5 C T A 3i ) P NC ( , 

T S °~ 2s c [X ± Y T sin (5c ±A 31 )+X T Y ± sin (5c T^ 31 )]±Y ± Y T sm2A 31 A ' 
= 5 _ [2X T s c - Y T cos (5 C ± A 31 )] P NC - [2X ± s c + Y± cos (5 C T A 31 )] P NC 

V± T ° 2s c [X ± Y T sin {So ± A 31 ) + X T Y ± sin (5 C T A 31 )] ± Sc; F ± r T sin 2A 31 1 ' 

It should be noticed that, given X± ~ C(l), Y± ~ C(e), and P^c ~ Pnc ~ C(e 2 ), C 
and T) are of order e 3//2 and e 1//2 , respectively. They are small compared to St (or sc) and 5 
which are of order e 1//2 and e°, respectively, justifying our perturbative treatment. Therefore, 
inclusion of the non-Cervera corrections does not produce sizable difference in the measured 
parameters. It is expected because of the smallness ~ e 2 of the non-Cervera corrections. 

Now let us discuss the parameter degeneracy. Since st — sc and 5t — 5c are small, the 
degeneracy solutions obtained with the Cervera et al. formula must give good approxima- 
tions to the ones obtained with our second order formulas. Therefore, we start from them; 
Suppose that all the degeneracy solutions (sa, 5a) (i = II, HI, ...VIII, reserving i = I for 
the true solution) are obtained by using the Cervera et al. formula. They are formally given 
by 

sa = fi(sci,5ci), 

5a = 9i(sci,5ci)- (31) 



The explicit expressions of the functions / and g are given in [24j]. Here, we take, as a 
concrete example, the case of the sign-Am 2 degeneracy whose solutions are labeled as III 
or IV. We mention how to extend our analysis to other types of degeneracies. Then, st — sc, 
the difference between si 3 obtained with our and the Cervera et al. formula, is given for the 
true / and the sign-Am 2 clone solution III as 

stiii — San — £,t( s ciii,5ciii) 

Stiii - 5 C iii = r] T (saii,5aii) (32) 
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Since sciii and Scm are given as functions of sci and 5c i, (|32|) with (l3Tj) parametrically 
solve st/// and <5rm as functions of sci and 5ci- For the intrinsic degeneracy solutions 
there is no need to flip the degeneracy sign in (1321) . For degeneracy solutions which involve 
flipping the octant of 6*23, X and Z in £ and 77 must be replaced by the ones with octant flip, 
X true -> X /aise = cot 2 9 23 X true and Z true ->■ Z /aZse = tan 2 9 23 Z true , as is done in Thus, 
all the degeneracy solutions obtained with use of our second order formula can be obtained 
by using the ones with the Cervera et al. formula as an intermediate step. The difference 
between the two is small for S13, order ~ e 3 / 2 , and somewhat larger, ~ e 1 / 2 for 5. 



VIII. INCLUDING THE NON-STANDARD INTERACTIONS 

The effects of possible nonstandard interactions (NSI) that would affect neutrino propa- 
gation in matter are usually described by adding the following additional term in the tilde 
basis Hamiltonian 



H 



NSI 



Ar A 



£ efj, 

e* e* i 



(33) 



The relationship between the i a p and th e e^ n (a, /3 = e,ft, r) parameters in the flavor basis 



is defined in (jSJ) 
of the order of e 



We assume, following [33 



the NSI elements £ a p (and hence e a p) are all 
It is a legitimate assumption because NSI comes from higher dimensional 
operators (dimension six or higher) which receives suppression of at least (Mw/Mnp) 2 ~ 
10~ 2 with new physics scale M^p- It should be mentioned, however, that the current bounds 
on the NSI parameters are quite loose, < 0.1 — 1 [66]. 

At small #i3 ~ e for which the e perturbation theory is applicable it was shown in 33 



that inclusion of NSI elements into the second-order formula can be done just by doing 
replacement (with slight change in the notations from (33| ) 



r&c 12 si2 -> r A c 12 s 12 + r A e efjL = E, 
s 13 e~ iS -> si 3 e" w + r A e eT = 6, 



(34) 



and nothing else, where e efl = (c 2 3,e eiJL — s 2 3,e eT ) and e eT = (s 2 3e efJi + c 2 3e er ). The z/ e -related 
oscillation probability is independent of e a/ 3 in the jJ, — t sector as well as e ee . Then, it is 
an interesting question to ask how the higher order corrections of Si 3 fit in into this picture. 
The second-order formula for P(v e — > v^) with NSI in our y/e perturbation theory can 
be obtained via a straightforward computation. The result is given by P(u e —> v^nsi — 
P{v e ->■ v^nsi-c + P{v e ->■ v^nsi-nc where 




sin 



(l-r A )AL 



(35) 
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+ 8s 



P{y e — » f^Nsi-Nc 



As 2 s 2 



13 , 



-I- 2s 2 S 2 
T ^*23*13 



1 - r A ) 2 
1 



'23 \^ee &tt 



l-r A ) 



2s 



ta 



13 



8C23S23S^ 3 |£ A(T | 



l-r A 



cos 



„2 (!+^) q„2 r A r A 
13 - S 



12 ^l-r^)3 



sm 



2 (l-r A )AL 



l-r A 



+ C 23 S23 £ur cos. 



sm 



;i-r A )AL 



1 - r A 



(AL) sin(l - r A ) AL 



AL\ r A AL (l-r A )AL 
sm sm , 



(36) 



where the NSI phase 4> a /3 is defined by e a p = \s a /3\^ al3 ■ Of course, 
obtained by the transformation C23 — >■ — S23 and S23 — > C23. 



v e — > v t )nsi 



can be 



The Cervera term 



contains the terms of order e 1 , e 3//2 , and e 2 . There exist two 



NSI dependent terms in P{y e 



Si3|£ e 



COS 



+ 5) or s 13 \e efl \ cos 



v^nsi-c m (l35|) of order e 3//2 which are proportional to 
Ai 1 Therefore, there is a confusion among the 



+ 2 



56, 671, but here 



three phases 5, (f) efl , and er which is known to exist in the small 9 regime 
with an amplified magnitude for large #13 ~ y/e. 

Notice that the NSI elements in the v^ — Vr sector as well as the diagonal ones are absent in 
the second-order formula of z/ e -related probabilities obtained with the e perturbation theory 
33J. The decoupling is supported for small $13 by the actual data analysis in which all 
the NSI elements are taken into account at the same time 56|. However, the formula ( 136]) 
tells us that it is no more true at order e 2 for large #13 of the order of the Chooz limit. It 
will require reconsideration of how to determine the NSI parameters simultaneously with 
the i/-mass enriched standard model parameters, the problem first addressed in [33[ whose 
elaboration is beyond the scope of this paper. 

Though the decoupling between the z/ e -related NSI elements and the ones in the — v T 
sector does not survive, a remnant still remains. The non-Cervera terms, P[y e — > v ^)nsi-nc 
in ( 136|) . consist only of order e 2 terms despite that the general theorems discussed in Sec. [TT] 
do not appear to apply to guarantee this property for the system with NSI. It is the 0(e 2 ) 
nature of the non-Cervera terms that leads to a remnant of decoupling, absence of ^-related 
NSI elements in the non-Cervera terms; 5 The terms induced by the substitution (I34p produce 
only higher-order terms of order ~ e 5//2 . 



IX. CONCLUSION 

In this paper, we have constructed a perturbative framework dubbed as "-y/e perturbation 
theory" to systematically compute higher-order corrections of S13 assuming that it is large, 
of the order of the Chooz limit. Despite a natural expectation that it could produce sizable 
corrections to the Cervera et al. formula, we have proven that they must be small, of 
the order of e 2 ~ 10~ 3 , where e = Am^/Am^. Nonetheless, we have observed that the 
correction terms nicely fill the gap between the exact oscillation probability and the Cervera 
et al. formula at baseline of several thousand km in a limited range of energy where they 
have enhancement due to the resonance effect. 



5 It can be argued that this property holds without doing any computations, as was done in the ArXiv 
version of this paper, arXiv:1103.4387vl [hep-ph]. 
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Possible large value of #13 may allow detection of the S dependent terms of ~ O (e 5 / 2 ) in 
future super-precision measurement because they are small only by a factor of ~ 5 compared 
to the O (e 2 ) terms. Therefore, we have computed the terms and found that they have 
different 6 dependence from the vacuum effect, in correlation to the matter dependent effects. 
However, it remains to be seen if these complicated and coexisting correlations can be 
resolved by actual experimental settings. 

Several characteristic features of the computed results prompted us to think about some 
general features of the 5 dependence of the oscillation probability. It resulted into the two 
general theorems stated and proved in Sec. [TT] and Appendix [A3 One of them (Theorem 
A) allows to understand why half-integral order terms in e are always accompanied with 
cosine or sine 5. While the other one (Theorem B) illuminates that the 5 dependence in 
the oscillation probability, both cosine and sine, must be suppressed at the small mixing 
angle or Am^/Am^ — >■ limits. Both of them cooperate to illuminate some characteristic 
features of the perturbatively computed oscillation probabilities 

We have investigated effects of the non-Cervera correction terms on parameter determina- 
tion and the degeneracy. Because the correction terms are small, their effect must be small. 
Nonetheless, the explicitly computed corrections to the degeneracy solutions obtained with 
use of the Cervera et al. formula may be of use when they are implemented into the analysis 



codes such as the ones described in 168, 69 



Finally, we gave a derivation of the second-order formula of the z/ e -related appearance 
probabilities with large #13 corrections in systems with NSI effects in propagation. The result 
of the Cervera terms is shown to be identical with the one obtained by the replacement to 
the generalized variables (1341) . producing enhanced NSI dependent terms of O (e 3 ^ 2 )- The 
decoupling of the NSI elements in the — v T sector from the ^-related probabilities is 
invalidated by the non-Cervera type corrections at order e 2 , which calls for reconsideration 
of the strategy for parameter determination. 



Appendix A: Theorem B for Suppression of CP Violation for Vanishing Mixing 
Angle 

Here, we attempt to prove the theorem B. We first describe our general strategy. We 
introduce a diagonal phase transformation matrix T = diag(e ia , e l/3 , e n ), and define a "hat 
basis" as v = Tv with the Hamiltonian H = THT*. We choose a, (3, and 7 such that H 
is free from any phases including 5. Since T is a diagonal phase transformation, it merely 
redefines phase of the flavor basis wave function, and hence it does not affect the probability. 
Therefore, if such choice of the phases is shown to be possible it implies that there is no 5 
dependence in the oscillation probability. 

We illustrate the proof by explicitly treating the case of #23 = case, because the situation 
is exactly the same in other cases. With the choice T = diag(l, 1, e~ lS ), the hat basis 



6 We thank Patrick Huber for communications for feasibility of implementing the analytic degeneracy 
solutions to achieve fast search for the degenerate minima, the possibility raised in j^ . 
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Hamiltonian can be written as H = H vac + dia.g(a/2E, 0, 0), where 

13 c 13 + ^31 c 13 s 13 

(Al) 



H V ac 



^2lSi2 c 13 + ^31 S 13 A21C12S12C13 — A 2 iS^ 2 Si3Ci3 + A 31 Ci3Si3 

A2lCi2S 12 Ci3 A 2 \c\ 2 — A21C12S12S13 

— A 2 iS^ 2 Si3Ci3 + A 31 Ci3Si3 — A21C12S12S13 A 2 iS^ 2 S^ 3 

+ A 3 ic? 3 



showing explicitly that 5-dependence disappear from the Hamiltonian by the T transforma- 
tion. Notice that here we have used the notation A^ = Am^/2E (ij = 21, 31), whose latter 
is different from the one in Sec. I VIII (We hope that no confusion arises.) Therefore, 5 goes 
away from the probability in the vanishing #23 limit. We can repeat the similar exercise for 
other vanishing limits of S12 with T matrix T = diag(e i<s , 1, 1). It is trivial to observe that 
there is no 5 dependence if #13 = 0. Therefore, we have shown that S dependence, both 
cos S and sin 5, disappears from the oscillation probability at the vanishing limit of one of 
the mixing angles for arbitrary matter density profiles. 

To really prove the theorem B we have to show that the 5-dependence goes away in the 
large mixing angle limit 9^ — > 7r/2 (ij = 12 and 23). One can show that the similar method 
works for c\ 2 —> and C23 —> with T = diag(e i<5 , 1, 1) and T = diag(l, e~ tS , 1), respectively. 
Assuming that the oscillation probabilities in all channels can be expanded into power series 
of Sij and the features stated above prove the theorem B for arbitrary matter density 
profiles. 

It is interesting to observe that the similar method fails for the limit C13 — > 0. It is 
perfectly consistent with the fact the factor C13 is missing in the coefficient of cos 5 term in 



the oscillation probabilities in the — v T sector (47|. However, since the suppression factor 
of the 5-dependent terms in the i/ e -related channels for constant matter density derived in 
the same reference |47| does contain c^ 3 , it is likely that the theorem B can be generalized 
to the one with suppression factor J instead of J r in these channels, a conjecture. 



Appendix B: S matrix elements to second order 

The S matrix elements can be computed by using (fT5|) and are written as sums over the 
terms of order ea, e 1 , e2, and e 2 . The T-conjugate elements can be obtained by Sp a (5) = 
S a p(—5). Notice that all the elements that are not shown below, except for the T-conjugate 
ones, are vanishing. 

= Sl3 e- iS —!— (e~ iAx - e~ irAAx ) (Bl) 

1 - Ta 




S£> = -c\ 2 r^iAx) 

§g> = sl 3 (3-^-) (iAx)e-^ - sisjY^ (e~ tAx ~ e+***) (B3) 
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The S matrix elements can be obtained from the S matrix elements by 



S 



S P 



q q 



q q q 

u [ie u fJ.iJ, u fj/r 



C23'S , ye + S23*SVe C^S^ + S 23 S TT + C23S23 {Sfj, T + S Tfl ) c 23*^Vr 



"23 ^A"" ~~ S 23^/1 + C23S23(5V r — S 
9 r> , 9 ri / rY , r> 



'23 

S 23^M + C 23 



23 ( -23' : '23V^rr ^/z/zj 

is^rr — C23S23(5'/xr + Sr/i) 



(B9) 



Appendix C: Second Order Expressions of Oscillation Probabilities 



In this AppendixOwe present some remaining expressions of the oscillation probabilities. 
In the v e —> v T channel, Per (i = 1, 3/2, 2) are given in each order by 



p(i) _ 4r 2 2 _ 
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